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FOREWORD 


This  investigation  was  assigned  as  a work  unit,  "Techniques  for 
Rapid  Subsurface  Exploration,"  of  an  RDT&E  Project,  "Permanent  Construc- 
tion Materials  and  Techniques  for  FY  70."  This  report  is  an  interim 
report;  a final  report  will  be  published  in  FY  74. 

Engineers  of  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  who  were  actively  engaged  in  the  investigation,  analysis,  and  re- 
port phases  of  this  study  were  Messrs.  R.  W.  Cunny,  Z.  B.  Fry,  R.  F. 
Ballard,  Jr.,  and  F.  K.  Chang,  of  the  Soils  and  Pavements  Laboratory 
(S&PL).  The  work  was  performed  during  the  period  December  1969  through 
March  1972  under  the  general  supervision  of  Mr.  J.  P.  Sale,  Chief,  S&PL. 
This  report,  which  is  essentially  the  same  as  a paper  presented  to  the 
42nd  Annual  International  Meeting  of  the  Society  of  Exploration  Geo- 
physicists, was  prepared  by  Messrs.  Ballard  and  Chang. 

COL  Levi  A.  Brown,  CE,  and  COL  Ernest  D.  Peixotto,  CE,  were 
Directors  of  WES  during  the  conduct  of  the  investigation  and  publication 
of  this  report.  Mr.  F.  R.  Brown  was  Technical  Director. 


CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 


British  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  units  as  follows : 


Multiply 

feet 

feet  per  second 
pounds  (mass ) 

pounds  (mass)  per  cubic  foot 
pounds  (force)  per  square  inch 


To  Obtain 


0.3048 

0.3048 


meters 


meters  per  second 


0.45359237  kilograms 


16.0185 


kilograms  per  cubic  meter 


0.6894757  newtons  per  square 
centimeter 
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SUMMARY 


The  objective  of  the  study  reported  herein  was  to  investigate  im- 
proved or  new  techniques  to  extend  the  depth  and  resolution  capability 
of  rapid  shallow-depth  (less  than  500  ft)  explorations  that  would 
quickly  and  economically  provide  an  engineer  with  accurate  information 
on  substrate  conditions.  This  report  does  not  advocate  the  adoption  of 
a new  test  technique  at  the  present  time . It  seeks  only  to  present  a 
discussion  of  potential  improvements  in  rapid  test  techniques. 

The  current  geophysical  techniques  for  rapidly  exploring  construc- 
tion sites  and  investigating  foundations,  i.e.,  electrical  resistivity, 
seismic  refraction,  and  vibratory  techniques,  are  discussed  herein. 

Preliminary  results  indicate  that  a study  of  surface  wave  phenom- 
ena can  measurably  extend  the  depth  limits  of  current  investigation 
techniques.  Utilizing  special  recording  and  interpretation  procedures, 
data  were  acquired  at  several  test  sites  and  compared  with  results  ac- 
quired by  more  conventional  methods . 

The  preliminary  results  of  a Rayleigh-wave  dispersion  technique 
have  yielded  velocities  consistently  greater  at  depths  of  50  ft  or  more 
than  the  conventional  vibratory  method.  These  apparent  differences  have 
not  yet  been  fully  resolved,  but  indications  are  that  these  discrep- 
ancies are  possibly  the  result  of  wave  distortions  caused  by  interfer- 
ing multiple  wave  trains  or  possibly  differences  in  stress  conditions  at 
the  source.  These  wave  patterns  and  resultant  group  velocities  will  be 
the  subject  of  future  study. 
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RAPID  SUBSURFACE  EXPLORATION 


REVIEW  OF  SELECTED  GEORfYSICAL  TECHNIQUES 

PART  I : INTRODUCTION 

Purpose  and  Scope 

1-  The  purposes  of  the  study  reported  herein  were  to  review  rapid 
subsurface  exploration  techniques  currently  employed  to  provide  engi- 
neering information  on  subsurface  geological  formations  that  may  be  used 
as  foundation  or  construction  materials,  and  to  investigate  alternate  or 
supplemental  techniques  as  improved  tools  for  rapid  subsurface  explora- 
tion. The  scope  of  this  report  involves  a cursory  review  of  existing 
rapid,  reliable,  low-cost  geophysical  methods  for  site  investigation  and 
location  of  construction  materials,  as  well  as  a detailed  description  of 
a new  Rayleigh-wave  dispersion  technique  for  improvement  of  existing 
rapid  exploration  capabilities . 

Background 

2.  The  strength,  compressibility,  and  permeability  properties  of 
soils  or  rocks  usually  control  the  design  of  foundations.  During  pre- 
liminary site  or  foundation  investigations,  rapid  means  of  indirectly 
determining  these  properties  or  indexes  to  the  properties  and  general 
stratification  over  large  areas  are  desired.  During  later  phases  of  the 
investigations  at  particular  sites,  more  detailed  or  specific  descrip- 
tions of  the  foundation  materials  are  desired.  The  subsurface  explora- 
tion techniques  used  to  determine  the  desired  properties  of  the  material 
can  be  divided  into  indirect  and  direct  techniques.1 

Indirect  techniques 

3.  In  geophysical  techniques  for  shallow-depth  exploration  (max- 
imum depth  to  bedrock  ~ 500  ft*),  the  seismic  refraction  and  the 

* A table  of  factors  for  converting  British  units  of  measurement  to 
metric  units  is  presented  on  page  vii. 
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electrical  resistivity  techniques  are  commonly  used.  ’ The  depth  as 
well  as  the  velocity  of  subsurface  materials  can  be  directly  determined 
by  a seismic  refraction  sounding;  the  dynamic  confined  compression  mod- 
ulus of  elasticity  can  be  calculated  from  the  velocity  if  the  density  of 
the  material  is  known  or  can  be  estimated.  The  electrical  resistivity 
technique  is  used  to  detect  the  variation  in  resistivity  of  earth  mate- 
rials, which  is  largely  dependent  upon  the  amount  and  salinity  of  the 
contained  water;  thus,  indirect  measurements  of  porosity,  saturation, 
and  permeability  are  possible.  The  seismic  refraction  and  resistivity 
techniques  each  have  certain  advantages  and  disadvantages  as  exploration 
tools.  Sometimes  the  two  techniques  complement  each  other  for  certain 
types  of  subsurface  investigations,  so  that  it  may  be  very  advantageous 
to  use  both  as  a means  of  achieving  speed,  economy,  and  reliability. 

Ten  years  ago,  the  U.  S.  Army  Engineer  Waterways  Experiment  Station 
(WES)  adopted  a combination  of  a vibratory  technique  in  conjunction  with 

seismic  refraction  technique  to  determine  the  in  situ  shear  and  Young's 
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elastic  moduli  of  soils.  The  vibratory  and  seismic  refraction  tech- 
niques were  used  to  measure  the  Rayleigh-  and  compression-wave  veloc- 
ities of  soils,  respectively.  The  gravitational  and  magnetic  tech- 
niques, generally  used  for  the  study  of  regional  geological  structure, 
are  without  depth  control  and  have  little  significant  engineering  value. 

4.  The  general  advantages  of  the  indirect  techniques  are: 

a.  Geophysical  techniques  of  exploration  identify  gross 
changes  in  character  of  subsurface  materials  by  surface 
measurement  of  the  changes  in  certain  physical  properties 
of  the  earth  at  or  near  its  surface . 

b.  Larger  areas  or  projects  of  greater  linear  extent  can  be 
explored  more  rapidly  and  economically  by  use  cf  geo- 
physical techniques  than  by  use  of  borings . 

c.  The  surface  measurements  as  determined  by  use  of  these 
techniques  indicate  average  conditions  within  a limited 
area  and  not  along  a single  vertical  or  inclined  line, 
such  as  is  the  case  when  the  boring  technique  is  used. 

d.  Data  acquired  by  use  of  these  techniques  can  be  used  to 
detect  the  subsurface  irregularities  that  often  are  missed 
by  borings . 


5.  The  disadvantages  are: 

a.  Although  the  general  character  of  subsurface  materials  can 
often  be  estimated,  the  materials  cannot  be  definitely 
identified  by  geophysical  techniques  alone;  therefore, 
these  techniques  must  be  supplemented  by  borings  in  which 
representative  and/or  undisturbed  samples  are  obtained. 

b.  The  interpretation  of  geophysical  data  requires  trained 
and  experienced  personnel. 

6.  Geophysical  techniques  are,  however,  especially  well  suited 
for  reconnaissance  and  preliminary  determination  of  for;aation  conditions 
of  large  areas  to  be  used  in  the  construction  of  missile  sites,  military 
underground  structures,  nuclear  reactors,  dams  and  reservoirs,  tunnels, 
highways,  airfields,  large  housing  projects,  etc.  Geophysical  tech- 
niques have  also  been  used  successfully  for  the  determination  of  water 
tables  and  deposits  of  gravel  or  other  engineering  construction 
materials . 

Direct  techniques 

7.  The  conventional  or  direct  techniques  utilized  for  subsurface 
exploration  and  determination  of  foundation  properties  include  boring; 
sampling;  measurements  of  penetration  resistance,  vane-shear  strength, 
plate-bearing  capacity,  and  California  Bearing  Ratio;  and  other  special 
purpose  tests  and  the  conduct  of  desired  laboratory  tests  on  disturbed 
or  undisturbed  samples  as  the  situation  requires.  The  direct  techniques 
of  exploration  are  relatively  slow  and  expensive  and  are  not  discussed 
further  in  this  paper.  The  authors  feel  that  the  thrust  of  this  paper 
should  be  directed  toward  a description  of  the  Rayleigh-wave  dispersion 
technique  as  a potential  new  tool  for  rapid  subsurface  exploration 
rather  than  discussing  "yesterday's  newspaper." 


PART  II:  RAYLEIGH-WAVE  DISPERSION  TECHNIQUE 


8.  It  is  generally  understood  that  the  conventional  seismic 
refraction  technique  will  provide  a compress ion-wave  velocity  and  an 
accurate  depth  determination  from  the  first  arrival  times  of  compres- 
sion waves,  and  that  the  shear -wave  velocity  at  depth  can  be  determined 
by  the  vibratory  technique,  which  has  a limited  depth  of  penetration  of 
about  150-200  ft.  Shear  modulus,  Poisson's  ratio,  and  Young's  modulus 
can  be  calculated  by  using  the  compression-  and  shear-wave  velocities 
and  the  mass  density  of  the  soil  or  rock. 

9.  Through  analyses  of  data  acquired  using  the  vibratory  tech- 
nique, it  has  been  found  that  the  waves  produced  by  the  vibrator  are 
predominantly  Rayleigh  waves.  It  is  considered  feasible  to  explore 
shallow  depth  by  means  of  Rayleigh  waves  generated  by  refraction  seismic 
tests  for  the  following  reasons : 

a.  Rayleigh  waves  as  well  as  compression  waves  are  generated 
by  the  refraction  seismic  tests. 

b.  Rayleigh  waves  are  more  easily  observed  than  shear  waves. 

c.  A pulse  traveling  with  Rayleigh-wave  velocity  should  be  a 
prominent  feature  on  any  seismic  refraction  record  unless 
the  instrumentation  has  been  designed  specifically  to 
reject  it. 

d.  Rayleigh  waves  generated  by  earthquakes  have  been  used 
successfully  as  a tool  for  exploration  of  the  earth's 
crust  and  mantle  structure  during  the  past  three  decades . 

10.  In  view  of  the  reasons  mentioned  above,  a proposed  seismic 
refraction  test  procedure  and  a plan  of  test  to  record  the  compression 
wave,  shear  wave,  and  Rayleigh  wave  on  the  same  seismogram  have  been 
considered.  Certain  preliminary  tests  have  been  conducted  for  compari- 
son of  Rayleigh-wave  velocities  determined  with  both  vibratory  and 
refraction  seismic  techniques . 

field  Procedure 

11.  The  field  procedure  for  recording  the  Rayleigh  or  R-waves 
during  refraction  seismic  tests  is  as  simple  as  that  for  recording  the 
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compression  or  P-waves.  Actually,  the  P-  and  R-waves  can  be  recorded  on 
the  same  seismogram  following  one  shot  if  the  instrumentation  incorpo- 
rates specially  designed  automatic  gain  control  (AUC)  circuitry.  Consid- 
ering that  a sharply  breaking,  large-amplitude,  first-arrival  P-wave  is 
necessary,  two  shots  must  be  made  with  less  sophisticated  instrumenta- 
tion: one  using  a large  amount  of  TNT  (3  to  5 lb)  for  recording  P-waves, 

and  one  using  a lesser  amount  of  TNT  (0.25  to  2 lb)  for  recording 
R-waves . The  size  of  the  charge  depends  on  the  shot-to-detector  dis- 
tance. At  a close  shot  distance  of  10  ft,  an  electrical  cap  may  be 
sufficient . 

12.  For  routine  dynamic  foundation  studies,  a twelve-geophone 
seismic  cable  with  25-ft  spacing  of  the  geophones  is  sufficient.  In  prac- 
tice, shots  are  fired  at  increasing  distances  from  each  end  of  the  cable 
spread  until  seismic  waves  from  the  bedrock  are  recorded  on  an  oscillo- 
graph. For  rapid  reconnaissance  in  a large  area,  a twelve -geophcne  seis- 
mic cable  with  geophones  at  50-ft  intervals  can  be  used.  For  rapid  sub- 
surface exploration,  two  seismic  profiles  of  fairly  long  length 
perpendicular  to  each  other  can  usually  be  completed  within  two  days . 

13*  The  instrumentation  system  employed  in  field  explorations  is 
similar  in  most  respects  to  that  used  in  shallow  refraction  explorations. 
The  band  of  frequency  response  including  the  geophone  amplifier  and 
galvanometer  is  wider  (around  the  range  of  2 to  200  Hz)  than  in  standard 
seismic  equipment.  The  geophones  are  of  the  velocity  type  and  are  espe- 
cially designed  for  low-frequency  measurement.  Present  WES  recording 
seismic  instruments  were  adapted  for  the  purpose  of  conducting  the  field 
experiments . 


Definitions  of  the 
group  and  phase  velocities 

14.  To  analyze  the  dispersion  characteristics  of  an  R-wave,  the 
group  or  phase  velocities  of  several  frequency  components  of  the  R-wave 
must  be  measured.  The  group  velocity  can  be  measured  by  assuming  that 
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the  wave  originates  as  a sharp  pulse  at  the  instant  of  the  shot  and  at 
the  surface  directly  above  the  shot  point.  On  this  assumption,  the 
group  velocity  v for  a travel  time  t of  a particular  cycle  of  pe- 
riod  T observed  at  the  minimum  horizontal  distance  X from  the  shot 
is  given  by 


(Reference  5)  (l) 


15.  In  the  case  of  soil  deposits,  the  seismic-wave  velocities 
increase  with  depth,  and  the  group  velocity  is  generally  less  than  the 
phase  velocity.  The  phase  velocity  can  be  found  by  numbering  the  peaks 
and  troughs  of  the  oscillations  and  following  them  frcm  one  adjacent 
trace  to  the  next  as  will  later  be  described  in  detail.  Assuming  that 
the  period  of  a given  cycle  in  the  sequence  does  not  change  drastically 
from  one  trace  to  the  next,  the  Rayleigh-phase  velocity  or  Rayleigh- 
wave  velocity  vr  can  be  approximated  by 


(Reference  5 ) (2 ) 


where  AX  is  the  distance  between  adjacent  seismometers  and  At  is  the 
apparent  time  required  for  the  R-wave  to  travel  this  distance.  (See 
fig.  la  for  graphic  representation  of  the  quantities.) 


Seismogram  Analyses 


l6.  Fig.  lb  is  a seismic  record  which  has  been  idealized  to  show 
the  prominent  wave  arrivals  at  two  adjacent  detector  locations  so  that 
the  R-wave  analysis  could  be  better  illustrated.  To  construct  R-wave 
dispersion  curves  of  phase  and  group  velocity,  the  following  steps  must 
be  undertaken : 

Phase  velocity 

a.  Number  the  crests  or  troughs  of  the  Rayleigh  wave  train 
for  all  traces . 

b.  Read  time  of  arrival  for  each  crest  (or  trough)  and  the 
average  period  T for  each  trace . 
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b.  IDEALIZED  SEISMIC  RECORD 

DEFINITIONS  UN  ORDER  OF  APPEARANCE  OR  USE): 

X = HORIZONTAL  DISTANCE  OF  SHOT  POINT  TO  FIRST  GEOPHONE  G1 
AX  = HORIZONTAL  DISTANCE  BETWEEN  GEOPHONES  G1  AND  G2 

t = RAYLEIGH-WAVE  ARRIVAL  TIME  (SUBSCRIPT  0 = FIRST  ARRIVAL, 
SUB-SUBSCRIPTS  1 AND  2 = GEOPHONE  NUMBERS) 

At  = APPARENT  TIME  REQUIRED  FOR  R-WAVE  TO  TRAVEL  Ax  DISTANCE 

0,1, 2. 3, 4 = SUCCESSIVE  TROUGHS  OF  RAYLEIGH  WAVE  TRAIN 

T = AVERAGE  PERIOD 


Vr  = RAYLEIGH-WAVE  VELOCITY  OR  V,  = 
L = WAVELENGTH  OR  L = VfT 

X 

Vg  = GROUP  VELOCITY  OR  vg  = J 


Ax 

At 


Fig.  1.  Example  of  seismogram  analyses 
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£.  Compute  the  incremental  time  At  between  adjacent  traces. 

d.  Determine  the  Rayleigh  phase  velocity  vr  for  each  wave 
number  by  vr  = AX/At  where  AX  = distance  between  ad- 
jacent seismometers. 

e.  Compute  the  wavelength  L for  each  wave  number  by 
L = vrT  . 

f.  Tabulate  data. 

£.  Plot  vr  as  a function  of  L or  approximate  depth,  L/2. 

Group  velocity 

a.  Number  the  crests  or  troughs  of  the  Rayleigh  wave  train 
for  at  least  one  trace . 

b.  Read  time  of  arrival  tQ  for  each  wave  number  and  its 
period. 

£.  Determine  group  velocity  Vg  for  each  wave  number 

- f (3) 

where  X - shot-to-detector  distance. 

d.  Plot  Vg  as  a function  of  T or  L (since  L = v T). 

Data  presentation  and  reduction 

17.  Providing  that  terrain  is  relatively  flat  and  background 
noise  levels  are  sufficiently  low,  the  oscillographic  data  presentation 
can  be  controlled  in  the  field  so  that  office  interpretation  time  will 
be  minimized.  This  can  be  accomplished  by  simply  adjusting  the  galva- 
nometer alignment  so  that  each  trace  will  be  spaced  in  proportion  to  the 
distance  between  adjacent  geophones.  In  this  way,  the  interpreter  may 
work  directly  from  the  field  records  if  he  so  desires.  The  field  data 
obtained  at  WES  test  site  2 (fig.  2)  can  be  used  to  illustrate  this 
and  several  other  interesting  points.  As  will  be  observed,  data  were 
recorded  in  four  225-ft,  12-trace  segments  from  geophones  that  were 
spaced  at  25-ft  intervals.  The  size  of  the  seismic  impulse  source 
ranged  from  electrical  caps  to  1/3  lb  of  C-4  explosive  charges.  (The 
size  of  the  charge  was  kept  to  an  effective  minimum  so  as  not  to  deface 
the  landscape.)  As  a means  of  positive  R-wave  identification,  traces  1, 
2,  and  3 were  the  registrations  of  a transverse,  radial,  and  vertical 
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Fig.  2.  Records  of  seismic  tests  at  WES  test  site  2 

geophone,  respectively,  located  at  the  first  station  of  each  of  the  four 
lines . 

18.  Utilizing  the  simultaneous  relationship  of  the  vertical  and 
horizontal  (radial)  geophones,  one  can  plot  the  trajectory  of  particle 
motion.  Fig.  3 shows  the  construction  of  these  diagrams  for  far-field 
stations  located  from  250  to  700  ft  from  the  seismic  source.  If  we  note 
that  both  upward  vertical  motion  and  outward  radial  motion  will  produce 
a downward  deflection  of  the  seismogram  trace,  one  should  expect  oscil- 
lations on  the  radial  trace  to  lead  those  on  the  vertical  trace  by  90 
degrees  in  phase.  This  is  seen  to  be  approximately  true  for  those  events 
labeled  Rayleigh  waves  on  the  record  shown  in  fig.  2.  Earlier  portions 
of  the  record  do  not  show  this  phase  shift  due  to  direct,  reflected,  and 
refracted  body  waves  causing  distortion.^  In  each  of  the  three  cases, 
the  particle  motion  roughly  takes  the  shape  of  a retrograde  elliptical 
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Fig.  3*  Trajectory  diagrams  of  particle 
motions  at  far-field  stations 

orbit,  thus  fulfilling  the  definition  requirements  for  classical  R-waves . 

19.  Observing  the  first  or  near-field  record  shown  in  fig.  2,  one 
can  readily  see  that  the  trace  signature  characteristics  are  completely 
different  from  the  far-field  records.  This  can  possibly  be  explained  by 
the  fact  that  the  distance  between  the  shot  and  the  geophones  was  less 
than  the  longest  wavelengths  that  were  observed  at  the  far-field  stations; 
therefore,  the  geophone  spread  is  too  close  to  observe  dispersion.  How- 
ever, a second  factor  must  also  be  considered;  the  seismic  source  for  the 
near-field  record  was  a single  electrical  cap  buried  only  0.5  ft  deep. 
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This  source  very  likely  possessed  insufficient  energy  to  excite  the  lower 
frequencies  associated  with  depth  penetration  and  consequently  greater 
soil  masses. 

20.  As  previously  stated,  interpretation  can  be  performed  directly 
on  the  records  if  so  desired.  For  example,  consider  the  second  record 

in  fig.  2 (250-475  ft)  for  illustrative  purposes.  Beginning  with  trace  3 
some  11  distinct  oscillations  can  be  observed,  numbered,  and  followed 
through  the  succeeding  traces  on  the  record.  A line  connecting  adjacent 
troughs  can  then  be  constructed  in  a manner  analogous  to  a conventional 
time-distance  plot.  Essentially,  the  reciprocal  slope  of  this  line  will 
yield  the  wave  velocity.  An  average  period  should  then  be  determined  so 
that  the  wavelength  can  be  computed  and  plotted  in  the  manner  previously 
described. 

Conversion  of  wave  veloc- 
ities  to  soil  properties 

21.  To  become  a useful  tool  in  terms  of  conventional  soil  prop- 
erties, the  P-  and  shear-  or  S-wave  velocities  determined  by  the  above 
method  must  be  related  to  elasti:  parameters  such  as  shear  modulus  G , 
Young's  modulus  E , and  Poisson's  ratio  v • This  conversion  can 
readily  be  accomplished  knowing  that  the  wave  velocity  is  dependent  on 
the  ratio  of  the  elasticity  of  the  medium  to  its  mass  density  p and  on 
the  wave  type . The  relationship  of  shear  modulus  G and  shear-wave 

velocity  v and  density  p is  as  follows: 
s 

G = v^p  (4) 

S-wave  velocity  and  surface  R-wave  velocity  are  related  by  Poisson's 
ratio.  For  a homogeneous  medium  and  for  Poisson's  ratios  ranging  from 
0.2  to  0.5,  the  difference  in  velocities  is  less  than  9 percent.  There- 
fore, for  practical  purposes,  S-waves  can  be  considered  to  have  the  same 
velocity  as  R-waves . 

22.  Poisson's  ratio  v of  a soil  can  be  determined  by  the  rela- 
tion of  S-wave  velocities  v and  P-wave  velocities  v 

s c 
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(5) 


G and  Poisson's  ratio 

E = 2(1  + v)G  (6) 

for  isotropic,  linearly  elastic  materials.  From  these  equations,  values 
of  E , G , and  v can  be  determined  through  the  measurement  of  S-  and 
P-wave  velocities,  provided  the  density  of  the  soil  is  known  or  can  be 
estimated. 

23.  Heukelom  and  Foster,  in  their  dynamic  testing  of  pavements 
using  the  vibratory  technique,  found  that  excellent  correlations  could 
be  estimated  with  known  depth  of  layers  if  they  empirically  assumed  a 
measured  velocity  applicable  at  a depth  equal  to  half  its  associated 
wavelength.  During  the  past  10  years,  WES  has  extensively  employed  this 
empirical  relationship  for  investigations  where  boring  data  were  also 
available  to  determine  depths  of  interfaces  for  different  materials  and 
agreement  has  generally  been  good;  consequently,  this  relationship  is 
assumed  valid  for  the  R-wave  dispersion  method. 

24.  Fig.  4 is  the  R-wave  velocity  for  WES  test  site  2 plotted  as 

a function  of  depth  (L/2)  prepared  using  the  R-wave  dispersion  technique, 
hi  addition,  data  obtained  by  using  the  surface  vibratory  technique  are 
shown  for  comparative  purposes.  It  is  interesting  to  note  that  excel- 
lent agreement  exists  between  the  vibratory  and  seismic  methods  t>  depths 
of  about  40  or  50  ft;  then  the  R-wave  dispersion  data  diverge  to  indi- 
cate higher  apparent  velocities.  This  can  possibly  be  explained  by  the 
fact  that  the  earliest  wave  arrivals  (generally  the  longest  wavelengths) 
are  a mixture  of  reflected  P-,  S-,  and  R-waves  that  tend  to  propagate  at 
a higher  apparent  velocity.  Unfortunately,  the  manual  method  of  deter- 
mining the  phase  velocity  is  seldom  accurate  enough,  and  it  is  generally 


Young's  modulus  E is  related  to  shear  modulus 
of  the  soil  by 
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Fig.  4.  R-wave  velocities  versus  depth  for  WES  test  site  2 

necessary  to  make  a Fourier  analysis  of  each  adjacent  pair  of  traces  and 
measure  At  from  the  phase  shift  of  the  Fourier  element  whose  period 

5 

is  T . Such  an  analysis  can  be  performed  with  available  computer  codes 
if  it  is  considered  necessary  by  the  data  user.  A second  possible  explan- 
ation is  the  difference  in  stress  levels  generated  by  the  two  methods  of 
testing.  Though  not  as  plausible  as  the  first  explanation,  it  must  be 
considered  as  a possibility  and  further  study  must  be  made  to  document 
possible  velocity  changes  as  a function  of  source  stress  level. 

25-  Once  the  elastic  moduli  have  been  determined,  they  may  also  be 
conveniently  displayed  as  a function  of  depth.  These  data,  which  were 
obtained  at  the  same  WES  test  site,  are  shown  in  fig.  5-  Once  these 

parameters  are  supplied  to  a knowledgeable  designer,  a structure  with 

7 

predictable  performance  can  be  designed.  Another  application  for  use 
of  these  soil  parameters  would  be  for  a dynamic  analysis  of  existing 
earthen  structures  such  as  dams  potentially  exposed  to  earthquake 
hazards . 
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Fig.  5*  Elastic  moduli  versus  depth  for  WES  test  site  2 
Theoretical  calculations 

26.  As  a means  of  supplementing  the  interpretation  of  data  ob- 
tained by  the  R-wave  dispersion  technique,  comparisons  with  theoretical 
dispersion  curves  based  on  soil  layering,  P-wave  velocity,  S-wave  veloc" 
ity,  and  density  of  the  media  can  be  made.  These  parameters  can  be 
determined  by  borehole  sampling,  conventional  seismic  refraction,  and 
some  engineering  estimates  of  material  properties.  Referring  again  to 
fig.  4,  one  can  see  the  comparison  of  theoretical  and  experimental 
values.  In  this  case,  all  three  methods  were  used  as  an  aid  in  estab- 
lishing layer  thickness  and  properties  of  the  WES  test  site.  The  model 
chosen  to  represent  this  site  was  representative  of  a four-layer  system. 
The  following  characteristics  were  assigned  tc  each  succeeding  layer: 


Depth 

Velocity 

J_fEs_ 

Wet  Density 

ft 

P-Wave 

S-Wave 

pcf 

0-12 

1,110 

350 

116 

13-56 

5,000 

650 

122 

57-76 

5,000 

980 

126 

77— 

10,000 

4550 

154 

Utilizing  these  assigned  descriptors  in  a R-wave  dispersion  computer 

g 

code  developed  by  Watson,  the  theoretical  fundamental  mode  for  this 
test,  which  was  treated  as  a four-layer  problem,  can  be  calculated.  It 
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is  apparent  from  fig.  4 that  some  differences  in  absolute  velocities  do 
exist . Future  research  will  seek  to  reconcile  these  differences . 


Evaluation 


27.  The  advantages  and  limitations  of  the  R-wave  dispersion  tech- 
nique are : 

Advantages 

a.  R-waves  are  in  most  instances  easily  observed  and  should 
be  prominent  features  on  any  seismic  refraction  record. 

b.  The  computation  of  phase  velocity  of  an  R-wave  is  easy. 

£.  The  R-wave  dispersion  technique  is  a combination  of  seis- 
mic refraction  and  vibratory  techi.^ques ; thus,  half  of 
the  field  working  time  could  be  saved  by  using  this  tech- 
nique. However,  it  must  be  stated  that  additional 
interpretation  time  spent  in  the  office  tends  to  minimize 
the  overall  cost  savings . 

Limitations 

a.  The  investigation  area  must  be  large  enough  that  adequate 
information  can  be  obtained.  For  example,  if  the  bedrock 
is  100  ft  deep,  the  length  of  the  profile  lire  (seismic- 
cable  spread)  should  not  be  less  than  500  ft. 

b.  Interpretation  of  the  R-wave  data  requires  the  services 
of  a competent,  trained  geophysicist. 

£•  The  R-wave  dispersion  technique  yields  velocities  that 

axe  consistently  higher  than  those  obtained  by  the  surface 
vibratory  technique  and  laboratory  test  methods. 


PART  III:  CONCLUSIONS 

28.  Preliminary  field  experiments  indicate  that  the  R-wave  dis- 
persion method  for  determining  shear -wave  velocity  by  surface  measure- 
ments is  indeed  feasible,  but  certain  questions  need  to  be  resolved. 

Most  apparent  is  the  lack  of  absolute  correlation  between  data  from  that 
method  and  surface  vibratory  data.  At  certain  test  sites  examined  at 
this  time,  the  R-wave  dispersion  method  yielded  consistently  higher 
velocities  at  depths  greater  than  50  ft.  Indications  are  that  the  dis- 
crepancy may  be  due  to  wave  interference  (model  superposition  on  the 
actual  seismograms)  or  possibly  stress  conditions  at  the  source.  How- 
ever, additional  study  is  needed  to  resolve  these  observed  differences. 
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